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The backreaction of nonlinear inhomogeneities to the cosmic expansion is analyzed in 
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the framework of general relativity with a cosmological constant. By defining the spatially 
averaged matter energy density, we find that the cosmological constant induces a new type 
of backreaction whose equation of state is p — — 4/3p, where p and p are the effective 
energy density and effective pressure of the backreaction in the averaged Friedmann universe. 
However, the effective density is negative, and thus it decreases the acceleration caused by 
the cosmological constant. 
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\ The observation of the isotropy of the cosmic microwave background radiation 

' (CMBR) 1 ) and large galaxy surveys, such as SDSS, 2 ) indicate that the universe is 

isotropic and homogeneous over scales of about 100 Mpc. However, local matter dis- 
tributions are highly inhomogeneous. This simple observation makes naive use of an 
isotropic and homogeneous model of the universe, namely the Friedmann-Lemaitre- 
Robertson- Walker (FLRW) model questionable, since the solution of the Einstein 
equations with an averaged homogeneous matter distribution does not solve the 
Einstein equations with a realistic matter distribution, because of its nonlinearity. 
Therefore, it is naturally expected that if we average a locally inhomogeneous uni- 
verse, the expansion of the averaged spacetime will be affected somehow by the local 
inhomogeneities. This effect is called the backreaction due to local inhomogeneities, 
and there have been many investigations of this backreaction. 3 )~ 8 )> n )~ 14 ) 

The problem of backreaction has attracted much attention recently. Observa- 
tions of Type la supernovae 15 )' 16 ) and the CMBR 17 )' 18 ) strongly suggest that the 
^ . cosmic expansion is accelerating. Understanding the source of this accelerated ex- 

pansion is one of the greatest unsolved problems in modern cosmology. 19 )' 22 ) Usually 
the existence of an unknown type of energy (dark energy, perhaps a cosmological 
constant) is assumed for the explanation of this acceleration. A possible alternative 
idea is that the energy caused by inhomogeneities leads to additional terms in the 
Friedmann equation, as if dark energy did exist. 23 )' 24 ) Although this possibility is 
unlikely, 25 )' 26 ) it is still important to derive a definite conclusion in regard to it and 
to carefully study its properties. 

In a previous study 10 ) we showed in the case without a cosmological constant that 
the nonlinear backreaction neither accelerates nor decelerates the cosmic expansion 
in a matter-dominated universe and that the backreaction behaves as a small positive 
curvature term. This conclusion was obtained by defining the averaged density of 
matter as a conserved quantity in the large comoving volume. Note that there are 
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many FLRW spacetimes that fit a particular inhomogeneous universe. Our choice 
is motivated by considering how one constructs a FLRW universe through actual 
observations, namely, by defining the evolution of the number density of galaxies by 
comparing with the comoving constant number density. It was also shown that the 
above results do not depend on the definition of the averaging procedure, and the 
results in the Newtonian gauge are consistent with those in the comoving synchronous 
gauge. 

However, the above results were obtained in the case with no cosmological con- 
stant, and it is not at all clear if the same conclusions apply in the case with a 
cosmological constant. On the basis of observations, it appears that the existence of 
the cosmological constant is very likely and thus it is worthwhile studying the back- 
reaction due to local inhomogeneities in a universe with a more general background. 
This is the purpose of the present study. 

We formulate the problem of the backreaction as generally as possible. For this 
purpose, we follow Kasai, Asada and Futamase 10 ) and employ the 3 + 1 formulation 
in general relativity. We restrict ourselves to the case of a vanishing shift vector, i.e., 
Ni = 0. The line element is written 

ds 2 = -{Ndtf + njdafda?. (1) 

The unit vector normal to a hypersurface foliated by t = const is denoted by 

^,0,0,o). (2) 

The extrinsic curvature is defined as 

= ^ tk %> (3) 

where the dot denotes the time derivative. The Einstein equations with a cosmolog- 
ical constant are reduced to the Hamiltonian constraint, the momentum constraint 
and (the trace of) the evolution equation as 

+ (K\f - JCjKU = 16irGE + 2 A, (4) 

K'^-Ki^^SnGJi, (5) 

K\ + NJCjK\ -N\\i = -AttGN(E + S) + NA, (6) 

where® -R is a 3-dimensional Ricci scalar curvature, ^ denotes the 3-dimensional 
covariant derivative, and the other symbols are defined as follows: 

E = T^nW = ^T 00 , (7) 
Ji = ~T vi n v = ^T 0l , (8) 

* r,r'>. (9) 
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We define the spatial 3-dimensional volume V of a compact domain D in a t = const 
hypersurface as 



V 

where 



/ Vld 3 x, (10) 
Jd 



7 = det( 7 y). (11) 

Here, V is considered to be a volume sufficiently large that we can assume periodic 
boundary conditions. 

The scale factor a(t) is defined in terms of the volume expansion rate of the 
universe: 

^ = v-vf D \f^ x - (12) 

Now we adopt the average procedure 

<A> = i jf A/yd 3 *. (13) 

Then, we have 

3^ = (NK\). (14) 
Furthermore, we define the deviation from uniform Hubble flow as 

V*; = NK^ - -J?,. (15) 

Then we can show (V l i) = 0. 

By averaging the Einstein equations, we obtain 

(£) 2 = 8 -^(N 2 E) - i(iV 2 MR) - \{{V\f - V)V\) + |(iV 2 ), (16) 

£ = -^f(N\E + S)) + i((V*i) a - V)V{) + ^(iViV 1 ; + NK\) + |<iV 2 >. (17) 

The geometrical treatment of these equations has studied by Buchert. 20 )' 21 ) 

In this paper we approximate the matter of the universe as irrotational dust 
whose energy-momentum tensor is given by 

T» v = pu^u u , (18) 

where is the 4-dimensional velocity of the fluid flow. 

Up to this point, our equations has been completely general, and a more detailed 
treatment requires explicit expressions for the inhomogeneous universe. In this pa- 
per we employ the post-Newtonian approximation for the metric to obtain a more 
concrete expression for the backreaction terms, which allows us to obtain a physical 
interpretation. 
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As the background solution, we employ the totally flat universe with a non- 
vanishing cosmological constant. It is then sufficient to consider the first-order cos- 
mological post-Newtonian metric, 

ds 2 = -(1 + 2(f)(x,t))dt 2 + a 2 (l - 2c/)(x,t))5 ij dx i dx j , (19) 

where 5ij denotes the Kronecker delta. Note that the potential (j> does depend on time 
in general. However, here we restrict ourselves to perturbations with scales much 
smaller than the horizon scale, Lh', otherwise the concept of averaging becomes 
unclear. This allows us to ignore the time dependence of the potential 4> in the 
Einstein equations. This approximation holds as long as i/Lu <C 1, where i is 
the scale of the perturbations we consider. The potential begins to decay when the 
cosmological constant begins to dominate the cosmic expansion. 
At first order, from the Einstein equations we obtain 

^ = 1^(24, + S)-^S, (20) 



8 ijV 1 = -2-(3a 2 -a 2 A)- 1 <l>, i , (21) 



where 5 = (p — Pb)/Pb is the density contrast, pf, is the background density, and 
contraction is taken with Sij. 

Next, we substitute these solutions into the right-hand sides of Eqs.(16) and 
(17), and retaining terms up to quadratic order in cf), the averaged Einstein equations 
become 

a\ 2 8ttG,_ , 1 , , , >A 



a) 3 W> + ^Mt) T <*>) 

a AttG , m 99 . 1 , , , . A . . 

- = 3-^00 + p b a 2 v 2 ) - stfWM + 3 . ( 23 ) 

where v 2 = 5ijV l vK We do not regard these equations as the Friedman equations for 
the averaged FRLW universe, and thus we do not use them to interpret the effect 
of the backreaction on the cosmic expansion. This is because these equations reveal 
that the background density pb is not conserved in the comoving volume, and thus 
pb is not appropriate as the mean energy density of the averaged FRLW universe. 

Now we define the conserved mean energy density for the averaged FRLW uni- 
verse. The mean density p is defined to be the quantity that satisfies the equation 

p = -3-p. 

a 

It is straightforward to derive the mean density using the energy conservation law 

Here, Eq.(21) has been used to rewrite the term containing (v 2 ) in terms of ((f>,i(f>,i), 
= Pb/Pcr, and f?A = A/3Hq, where Hq and p cr = Hq/8ttG are the Hubble 
parameter and critical density, respectively. 
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We rewrite Eqs. (|22l) and (123j) in terms of the mean density p: 

"P+T-n 77 + 77^ <^,i)» (25) 



a / 3 



* = ™ f+ A-!jgL M< i>. (26) 

a 3 3 bJ/ m 

These are the Friedmann equations for the averaged FLRW model of a locally in- 
homogeneous universe. The third terms on the R.H.S. of the above equations are 
interpreted as the non-linear backreaction due to local inhomogeneities. There are 
two types of backreaction. One is proportional to a~ 2 , and hence behaves as a 
small positive curvature term in the Friedmann equation. This is the same as in the 
Einstein-de Sitter case. 10 ) The other is proportional to a, and hence behaves as a 
phantom fluid in this sense. However, this term acts to decelerate the cosmic expan- 
sion, as shown in Eq.(26). The backreaction induced by the cosmological constant 

has strange properties. If we define an effective energy density and effective 

pressure P^ of the backreaction, then we find 

^ - A<*<** < 28 » 

Thus, the equation of state is P^ = — . 

When the cosmological constant begins to dominate the expansion, the potential 
begins to decay. Thus the backreaction also decays. Note that we have ignored the 
time dependence of the potential when averaging the Einstein equation. As discussed 
above, this is allowed as long as we consider perturbations with scales much smaller 
than the horizon scale. In fact, we can evaluate the backreaction using the power 
spectrum P^(k) of the gravitational potential as 

''' k k 2 P<i>{k). (29) 



(2k) 



Employing a realistic ylCDM model for the power spectrum, 27 ) we have confirmed 
that changing the minimum momentum in the integration from i^o to OAHq results 
in a change in the value of {(j) t i<fi,i) that is less than 0.3%. 

We have found that the gravitational interaction between local inhomogeneities 
and the cosmological constant induces a peculiar backreaction which appears in or 
near the epoch of the structure formation only. Its effects on the global expansion and 
the distance-redshift relation are in general very small. For example, the backreaction 
changes the luminosity distance from the present to z = 1000 by less than 0.001% 
from the standard distance in the /1CDM model with fi m = 0.3 and Qa = 0.7. 
However, it is interesting that the effect of the backreaction changes the observational 
effect if we properly take into account the large scale perturbations which we have 
ignored in this letter. 
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Also, the local inhomogeneities are important in cosmological observations. In 
fact, the scattering of light rays from distant sources by local inhomogeneities causes 
a redshift-dependent dispersion in the distance-redshift relation which will be essen- 
tial to properly interpret cosmological observations such as the magnitude-redshift 
relation of Type la supernpovae 9 ) and the evolution of the quasar luminosity func- 
tion. 
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